Quasi-discrete microwave solitons in a split ring resonator-based left-handed coplanar waveguide 
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We study the propagation of quasi-discrete microwave solitons in a nonlinear left-handed coplanar waveguide 
coupled with split ring resonators. By considering the relevant transmission line analogue, we derive a nonlinear 
lattice model which is studied analytically by means of a quasi-discrete approximation. We derive a nonlinear 
Schrodinger equation, and find that the system supports bright envelope soliton solutions in a relatively wide 
subinterval of the left-handed frequency band. We perform systematic numerical simulations, in the framework 
of the nonlinear lattice model, to study the propagation properties of the quasi-discrete microwave solitons. 
Our numerical findings are in good agreement with the analytical predictions, and suggest that the predicted 
structures are quite robust and may be observed in experiments. 



PACS numbers: 41.20.Jb, 42.65.Tg, 78.20.Ci 

I. INTRODUCTION 

Artificially engineered metamate rials have unique electro- 
magnetic (EM) properties, which are of great interest both 
from the basic physics viewpoint and for numerous applica- 
tions 0J[3]]. In such media, the effective permittivity e and per- 
meability /i is such that, in certain frequency bands, the meta- 
material displays either a right-handed (RH) behavior (e > 0, 
fi > 0) or a left-handed (LH) behavior (e < 0, fi < 0); in 
other words, energy and wavefronts may travel in the same 
or opposite directions in the RH or the LH frequency re- 
gion, respectively. Metamaterials with a LH behavior, namely 
LH-metamaterials (LHM), exhibit negative refraction at mi- 
crowave HEl or optical frequencies 0. 

Apart from the effective medium description, there has also 
been a large interest in studying equivalent transmission line 
(TL) models for LHMs in the microwave frequency region. In 
such a case, the effective permittivity e and permeability \i are 
directly connected to the serial and shunt impedance of the 
TL model; this way, so-called, composite right-left handed 
(CRLH) TLs (7 | may exhibit either a RH or a LH behavior 
depending on the frequency band. In practice, CRLH-TLs 
may be implemented with the coplanar waveguide (CPW) [8 1 
or microstrip technology [9|. Importantly, CRLH-TLs have 
led to many microwave applications and devices, including 
dual-band branch-line couplers, asymmetric backward-wave 
directional couplers, resonators, antennas, and so on UKUEJ. 

On the other hand, nonlinear metamaterials, namely struc- 
tures in which e and /i (or the serial and shunt impedance in 
the respective TL models) depend on the EM field intensi- 
ties (or voltages and currents in the TL models), have also 
been a subject of interest. Such structures may be imple- 
mented by embedding an array of wires and split ring res- 
onators (SRRs) into a nonlinear dielectric iflOl fTTI . or by in- 
serting diodes into resonant conductive elements (such as the 
SRRs) |[T2lfl5l . Nonlinear metamaterials may prove useful 
in various applications, including "switching" the material 
properties from left- to right-handed and back, tunable struc- 
tures with intensity-controlled transmission, negative refrac- 
tion photonic crystals, etc. Furthermore, fundamental effects, 



such as harmonic generation, nonlinearity-induced localiza- 
tion of EM waves and soliton formation, are possible. More 
specifically, relevant nonlinear phenomena [[TBI [TTl and soli- 
ton formation have already been predicted to occur in nonlin- 
ear metamaterials, using either the effective medium descrip- 
tion (see, e.g., Refs. [ 18-25]) or the TL description (see, e.g., 
Refs. 1 26 -28]). From the viewpoint of nonlinear TL experi- 
ments, pulse propagation |29| and envelope soliton formation 
Poll were recently observed (see also the review of Ref. PTTl 
and the relevant work of Ref. |32|); moreover, analytical ap- 
proximations, based on a continuum nonlinear Schrodinger 
(NLS) equation, allowed the description of bright l32l [331 or 
dark [32, 34] envelope solitons observed in the experiments. 

In this work, we study quasi-discrete microwave solitons 
that may be formed in a planar LHM. Such structures have 
the advantage of being easily fabricated (by means of stan- 
dard mask/photoetching techniques), they are compact, and 
are compatible with monolithic microwave integrated circuits. 
At this point we should mention that, generally, metamaterial 
TLs, operating at microwave frequencies, are artificial lines 
consisting of a host line loaded with reactive elements. Such 
TLs can be implemented by means of two main approaches: 
(i) the so-called "CL-loaded approach", where RH-TLs are 
loaded with series capacitances and shunt inductances [7-9] 
(see also Ref. Il32l for results in a nonlinear case), and (ii) the 
so-called "resonant-type approach", where the TLs are loaded 
with sub-wavelength resonators, such as SRRs [35, 36 1. TLs 
of the latter type exhibit controllable electrical characteris- 
tics, beyond what can be achieved in conventional TLs, im- 
plemented, e.g., in printed circuit boards (PCBs): the size of 
such a TL is determined by the size of the resonators and, thus, 
they can easily be miniaturized. Note that the fabrication of a 
pertinent prototype device, based on a CPW with an array of 
SRRs being etched at the bottom of the substrate, was first in- 
troduced in Ref. [35 1, while the corresponding TL model was 
presented in Ref. Il37ll . 

Here, we consider a nonlinear counterpart of the TL model 
of Ref. 11371 . with the shunt capacitors of the model being 
nonlinear; such a nonlinear TL model may be implemented 
by incorporating a nonlinear dielectric thin film in the struc- 
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FIG. 2: The unit cell circuit of the SRR-CPW model. 



ters of the SRR-CPW model, the second analyzes the method 
of multiple scales used to obtain the relevant NLS equation, 
while the third touches upon the continuum limit of the quasi- 
discrete approximation developed herein. 



FIG. 1: (Color online) Top panel: Configuration of the SRR-CPW 
structure: the dielectric layer of relative dielectric constant e r \ and 
width hi — hi is linear, while the one of relative dielectric constant 
e r 2 and width hi is nonlinear. Bottom panel: Layout of the bottom 
plate of the configuration, where the SRRs are placed. Shown are all 
relevant dimensions appearing in the equations of Appendix A. 



ture, whose dielectric constant may be controlled by a proper 
bias voltage (see, e.g., Refs. [38-40|). In our analysis, start- 
ing from the discrete lump element model of the CRLH-TL 
under consideration, we derive a nonlinear lattice equation, 
which is then treated in the framework of the quasi-discrete 
(alias quasi- continuum) approximation (see, e.g., Ref. ATI 
and Ref. Il42l for a review): this way, seeking for envelope 
soliton solutions of the nonlinear lattice model, characterized 
by a discrete carrier and a continuum envelope, we employ 
a multi-scale perturbation method to derive an effective NLS 
equation. The coefficients of this equation, which determine 
the type (bright or dark) of the envelope soliton, are found and 
it is shown that bright NLS solitons are supported by the SRR- 
CPW structure in a relatively wide range of frequencies inside 
the LH frequency band. Our analytical predictions are corrob- 
orated by numerical simulations, which reveal (apart from the 
basic properties) the robustness of the predicted quasi-discrete 
microwave solitons. 

The paper is organized as follows: in Sec. II we present the 
SRR-CPW nonlinear structure and derive a nonlinear lattice 
model describing the evolution of the voltage across this struc- 
ture; then we employ the quasi-discrete approximation and de- 
rive the NLS equation describing envelope solitons that can be 
supported in this setting. In Sec. Ill we present numerical sim- 
ulations in the framework of the nonlinear lattice model, re- 
veal the propagation properties of the predicted quasi-discrete 
solitons, and compare the numerical results with the analyti- 
cal predictions. We also propose changes in the configuration 
to increase the domains of existence of solitons. Finally, in 
Sec. IV we discuss our conclusions. Our presentation is com- 
plemented by three appendices: the first details the parame- 



II. THE MODEL AND ITS ANALYTICAL 
CONSIDERATION 

A. The nonlinear lattice model 

We start by considering a nonlinear version of the SRR- 
CPW model introduced in Refs. 0511371 . as shown in the top 
panel of Fig. |2| here, the structure incorporates a nonlinear 
dielectric film, of relative dielectric constant e r2 and width 
h,2- The rest of the SRR-CPW configuration is identical to the 
one of Refs. 11351 [371 : one may observe the SRRs (of exter- 
nal radius R) at the bottom plate of the structure, which are 
aligned with the slots (of width G and separation distance D) 
at the top plate of the structure. Notice that below we will 
focus on the case where the nonlinear dielectric film is intro- 
duced as shown in the top panel of Fig. |2j according to this 
consideration, nonlinearity is only introduced in certain ele- 
ments (i.e., shunt capacitors) in the equivalent discrete unit- 
cell model of the system. Nevertheless, below we will first 
consider the case where the serial capacitors, associated with 
the SRRs, are also nonlinear; in practice, this can be done, 
e.g., by inserting diodes in the SRR slots lfT2lfl5l . 

The discrete element model (unit cell) of the considered 
SRR-CPW structure is shown in Fig. [2] Here, Lr and Cr 
denote the equivalent per section inductance and capacitance 
of the line, respectively, L s and C s are the equivalent induc- 
tance and capacitance of the SRR, which is coupled with the 
transmission line, while the inductance is the equivalent 
inductance of the shunt strips. The above elements are di- 
rectly connected with the physical parameters of the SRR- 
CPW structure (see, e.g., Ref. [43] and details in Appendix 
A). 

Let us now consider Kirchhoff's voltage and current laws 
for the SRR (see the L S C S combination in Fig. [2]) equivalent 
circuit, namely, 



di 2 dJ n 

L 's— J T — Vn-1 ~ v n — L R ——, 

at at 

d ( dJ n 
11 = s dt 1 n_1 ~ Vn ~ R ^dJ 



(1) 

(2) 



The above equations, together with the auxiliary Kirchhoff's 
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current law 11+^2 = J n , lead to the following system for the 
unknown voltages V n and U n : 



cies related to the SRR and the RH part of the circuit, respec- 
tively; the above frequencies read: 



L R ^[C s (U n )U n -C s (U n+1 )U n+1 ] fse = 

+ (l + j^) (U n - U n+1 ) + (V n -i - 2V n + V n+1 ) = 0, fa = 

(3) 
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2tt^/L s C s 



2ttVLrC q 



(10) 



L R —[C(V n )V n ] - (K-i - 2V n + V n+1 ) 

+^V n -{U n -U n+1 )=0, (4) 

where C = 2C R and L = Li/2. At this point, in order 
to further simplify Eqs. Q-Q, we need to make some con- 
siderations concerning the nonlinear (voltage-dependent) ca- 
pacitance C = C(V n ) and C s (U n ). As discussed above, C 
is effectively filled by a nonlinear dielectric, whose dielectric 
constant may be controlled by a proper bias voltage, say Vb 
(see, e.g., Refs. [38-40 1 and discussion below). On the other 
hand, the effective nonlinearity of C s is due to the presence 
of a diode, biased at a constant voltage, say Ub- Assuming 
that the voltages V n and U n , respectively applied in the ef- 
fective capacitors C and C s , do not change significantly from 
their relevant bias voltages Vq and Uq, we can Taylor expand 
C(V n ) and C s (U n ) taking into regard only the lowest order 
terms, namely, 

C(V B + V n ) « Co + C' V n , (5) 

C s (U B +Un)^C s o+C' s0 U n , (6) 

where Cq = C(Vb) and C s q = C s (Ub) are constant capac- 
itances corresponding to Vb and Ub, respectively, and sim- 
ilarly C = {dC/dV) \ Vb , C s0 = (dC s /dU) \ Ub , take 
constant values depending on the particular forms of C(V) 
and C S (U) (see also below). Next, substituting Eqs. (|5j— 
(|6]l into Eqs. Q-Q, and introducing the scale transforma- 
tions t — > uj s ht [where oj 2 h — (LCo) -1 ] and {V n ,U n } — > 
[C (2Co) _1 ]{14, U n }, we cast the system of Eqs. into 
the following form: 



d 2 V n 



dt 2 

= (3 2 (U n - r, 



P 2 {V n+1 + Vn-l - 2V n ) + V n + 



d'Vn 2 

dt 2 



n+1), 



(7) 



and 



(Un ~ U n+1 ) + J 2 (U n ~ U n+1 ) + °^ 2 (Ul ~ Uf l+1 ) 



dt 2 

+( 7 2 -m 2 )(K+i + K- 



2V n ) = 0, 



where the constant parameters are given by: 

CqC s0 f se fr/.H 

°'=7v 7 ^, 7=7-, P=-f — > f 1 = 
L'O^sO Jsh Jsh 



fsh 



(8) 



(9) 



In the above expressions, f se and f^ denote series and shunt 
frequencies, while f s and fnH denote characteristic frequen- 



The system of Eqs. |7]i-([8]l is a model of two coupled non- 
linear lattice equations, for the unknown voltages V n and U n , 
describing the dynamics of the system. The analysis of the 
above system by means of analytical (and/or even numerical) 
techniques is a far more involved task and will be deferred 
to a future study. In fact, we choose to consider a variant 
of the above model, which is analytically (and numerically) 
more tractable, corresponding to the case where only the ca- 
pacitance C (associated with the CPW structure) is nonlin- 
ear, while the capacitance C s is linear (i.e., nonlinearity in the 
considered SRR-CPW structure is only introduced by the in- 
sertion of the nonlinear dielectric film, as shown in Fig.[T]i. In 
such a case, the parameter o (governing the SRR nonlinear- 
ity) becomes a = 0, and the system of Eqs. <|Vj-<[8j can be 
decoupled, leading to a single nonlinear lattice equation for 
the function V n , namely: 



^Vn , ,„ , 2 ,d 2 V n 



+ (1 + 7^)- 



+ r 



dt 4 

M 2 /3 2 (K+i + V n -i - 2V n ) + i 2 V 
,d 2 V 2 



dt 2 



^=0 

dt 4 



2V n ) 



(11) 



It is now useful, for the purposes of our analytical and nu- 
merical considerations, to adopt experimentally relevant pa- 
rameter values. First, following Ref. Il37l . the parameters re- 
lated to the CPW structure are chosen as follows. The waveg- 
uide has a width a = 23.7 mm and a thickness hi = 1.27 mm, 
the central strip width is D = 7 mm, the width of the slots is 
G = 1.35 mm, while the main - linear - dielectric substrate 
(namely a Rogers RO3010), of width hi - h 2 = 1.268 mm, 
has a relative dielectric constant e rl = 10.2. The values for 
the SRR characteristics are also borrowed from Ref. Il37l and 
are assumed to take the following values. The internal radius 
is r = 2.4 mm, the distance between the rings is d = 0.2 mm 
and the rings width is w = 0.6 mm. As far as the nonlinear 
dielectric film is concerned, following Ref. Il40l . we have as- 
sumed a strontium barium titanate (SBTO) paraelectric thin 
film, of width h 2 — 2 /im, and relative dielectric constant 
e, 2 = 300. Using the above parameter values, we may deter- 
mine (as per the relevant equations provided in Appendix A) 
the values of the effective capacitances and inductances in- 
volved in the SRR-based CPW structure. This way, we find 
that L R = 4.11 nH, L = 0.9 nH, L s = 1.33 nH, C s = 4.9 pF, 
and Co = 2.44 pF. Notice that the value of Co is also obtained 
consistently from the effective voltage-dependent capacitance 
of Ref. 1 40 1, namely: 



C(V) = C 



1 



(bo + bM+hV 



(12) 
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where C Q = 1.5 pF, b = 0.49, bi = 0.25 V" 1 , while V B = 
4.38 V is the constant (DC) bias voltage. According to the 
above expression, the values of the constant coefficients in 
Eq. ^ are given by C = C Q [1 + [6 + 2b 1 V B )- 1 } = 2.44 pF 



[identical to the result of Eq. (Al i], and C = C b\{h 



2hV B y 2 = 0.147 pF/V. 

According to the above, the frequencies in Eq. (IOi, 
which are in the microwave regime, take the values f se = 
2.272 GHz, f sh = 3.395 GHz, f s = 1.975 GHz, and 
Jrh = 1-592 GHz; accordingly, the values of the normalized 



parameters 7, /i and /3 appearing in Eq. (Hi take the following 
values: 



7 = 0.66, n = 0.58, (3 = 0.47. 



(13) 



Thus, in our simulations (see Sec. III. A), we will use the above 
values to investigate quasi-discrete microwave solitons in this 
setting, and also discuss possible modifications of the consid- 
ered setup in order to study how relevant changes in param- 
eter values affect the domains of existence of these nonlinear 
structures (see Sec. III.B). Notice that for the above mentioned 
choice of the physical parameters, the time unit associated to 
Eq. ( 1 1 1 is to = (2vr/ s /j) _1 « 50 picoseconds, while the volt- 
age unit is vo = 2Cq/C « 33 Volts. 

Before proceeding further, we should note the following. 
As indicated by Eq. |5]), in our considerations we take into 
regard only the first order approximation in the C(V n ) depen- 
dence, while this restriction may not be accurate enough for 
artificial nonlinearities (induced, e.g., by inserting nonlinear 
elements, such as diodes, in the SRR slots - see Ref. |[T3l 
and ll32l for relevant theoretical and experimental studies, re- 
spectively). Nevertheless, in our case, where the nonlinear- 
ity is induced by the insertion of the nonlinear dielectric thin 
film (see Fig. [TJ, the considered approximation is quite rea- 
sonable: indeed, in our numerical simulations (see Sec. Ill) 
we use a value for the initial voltage equal to Vq = 0.5 V 
(in physical units). For such a value of the voltage (simi- 
lar, and even smaller, values have also been used in relevant 
experimental works [29-34)), it can be found that C(Vq) = 
C{V B ) + C' Q V = (2.44 + 0.0735) pF = 2.51 pF, i.e., « 3% 
higher than the value of Co = 2.44 pF. If we had taken into 
regard the quadratic term in the Taylor expansion in Eq. d3), 
namely ^C V 2 , this term would take the value 0.00575 pF, 
i.e., only s» 0.2% higher than the above mentioned value of 
C(Vb) = 2.51 pF (that corresponds to the lowest order of 
approximation). Thus, according to the above arguments, and 
given that the initial voltage value Vq can also be controlled by 
other additional parameters stemming from our analysis [see 
parameters e and r\ in Eq. ( |23j ) below], we will proceed by 
analyzing Eq. (Hi that takes into account only the first order 
approximation in the C(V n ) dependence [as perEq. |5}]. 



B. The quasi-discrete approximation and the NLS model 

In this Section, we will employ the quasi-discrete approx- 
imation (see, e.g., Refs. fiTl l42l and Appendix B). Gener- 
ally, this approach is a variant of the multi-scale perturba- 
tion method, which is a well-known powerful tool to derive 



effective evolution equations (valid under certain conditions 
and for appropriate spatial and temporal scales) that are much 
simpler than the original models |45ll . In our case, since our 
original model [cf. Eq. (Ill] is actually a nonlinear dynami- 



cal lattice, we adopt the quasi-discrete approximation due to 
the fact that it takes into regard the discreteness of the system: 
this approach, allows for the description of quasi-discrete en- 
velope solitons (satisfying an effective NLS model), which 
are characterized by a discrete carrier and a slowly-varying 
continuum pulse-like envelope. Notice that, alternatively, one 
could adopt a continuum approximation, i.e., take the contin- 
uum limit of Eq. (Ill and analyze the latter in the framework 



of a multi-scale perturbation scheme (as in Refs. Il27l[33l[34ll ). 
However, as we will show below (and as was the case in 
Ref. l32l ). the quasi-discrete approximation is more accurate 
than the continuum one on providing estimates for the do- 
mains of existence of envelope solitons. Furthermore, the an- 
alytically determined soliton profile and characteristics (such 
as the center of mass and width) will be found to be in good 
agreement with direct simulations obtained in the framework 
of the original lattice model of Eq. ([TT| - see Sec. III. A below. 
We seek for solutions of Eq. \\ \\ in the form: 



(14) 



where Vt (£ = 1, 2, • • •) are unknown envelope functions de- 
pending on the slow scales X = e(n — v g t) (where v g is the 
group velocity, to be determined in a self-consistent manner) 
and T = e 2 t; here, < e 1 is a formal small parame- 
ter related to the soliton amplitude (see below). Additionally, 
the function exp(i9 n ), with n = cut — kn (with u> and k 
denoting frequency and wavenumber, respectively) describes 
the carrier. In the above ansatz, the envelope (carrier) is ob- 
viously continuous (discrete) in space; the results obtained in 
the framework of the quasi-discrete approximation, may be 
directly viewed in the continuum limit (and would correspond 
to the continuum approximation) of k — > 0; see also Appendix 
C. 



Substituting Eq. ( 14 1 into Eq. ( 11 1 we obtain the following 
results (see more details in Appendix B). First, to order 0(e) 
(linear limit), we derive the following dispersion relation: 

w 4 -(l+7 2 +4/3 2 sin 2 ^V 2 +4/3Vsin 2 ^+ 7 2 = 0. (15) 



In Fig.|3]we plot the frequency / = Lo/2i:f s h (in GHz) as 
a function of the wavenumber k (in rad/cell). It is clear that, 
apart from the shown gaps (where EM wave propagation is not 
supported), there exist two different frequency bands (where 
propagation of EM waves is possible): a high-frequency (HF) 
and a low-frequency (LF) band, with 3.395 GHz < / < 
4.734 GHz and 2.102 GHz < / < 2.272 GHz, respectively. 
It should be noticed that the lower (upper) cutoff frequency 
of the HF (LF) frequency band is equal to f s h (f se ) - see 
Eq. (10 1. Obviously, in the HF (LF) band the product of 
the wavenumber k and the group velocity v g = dui/dk is 
positive (negative) and, thus, energy and wavefronts travel 
in the same (opposite) directions in the HF (LF) frequency 
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0.5 1 1.5 2 2.5 3 
Wavenumber (rad/cell) 

FIG. 3: (Color online) The dispersion relation, showing the fre- 
quency / (in GHz) as a function of the wavenumber k (in rad/cell). 
Solid (blue) lines and dashed (red) lines correspond to the quasi- 
discrete and continuum approximations [cf. Eqs. {T3} and «[CTJ, re- 
spectively]. There exist a high-frequency RH band, for 3.395 GHz < 
/ < 4.734 GHz, and a low-frequency LH band, for 2.102 GHz < 
/ < 2.272 GHz; the inset shows a magnification of the dispersion 
relation in the LH frequency band. 



region. Thus, in the LF band (which is clearly a LH fre- 
quency band) the considered SRR-CPW structure apparently 
behaves as a LH transmission line; the dispersion relation 
in this LH frequency band, is shown in the inset of Fig. [3] 
Note that in the continuum limit (see pertinent results, for 
k — » 0, in Appendix C) the LH frequency band under con- 
sideration becomes slightly wider, i.e., it extends in the inter- 
val 1.976 GHz < / < 2.272 GHz; a similar effect is also 
observed for the RH band (/ > 3.395 GHz). 

Next, proceeding to the next order [0(e 2 )] in the perturba- 
tion scheme, we obtain the group velocity, given by: 



duj 
~dk 



p 2 smk{uj 2 -^ 2 ) 



2oj 3 - (1 + 7 2 +4/3 2 sin 2 



(16) 



Finally, to order 0(e 3 ), we obtain a nonlinear evolution equa- 
tion for the unknown voltage V\ (X, T), namely the following 
NLS equation, 

%d T Vi + Pd\V x + Q|Fi| Vi = 0, (17) 

with dispersion and nonlinearity coefficients, P and Q respec- 
tively, given by the following expressions: 



P 



d 2 UJ 

dk 2 



(1 + 7 2 + 4/3 2 sin 2 



6lo 2 )v 2 



-4w 3 + 2(1 + 7 2 + 4/3 2 sin 2 |)w 



4/? 2 w s o;sinfc + /3 2 (lj 2 — [i 2 ) cosfc 
-4w 3 + 2(1 + 7 2 + 4/3 2 sin : 



2 fcN 

2 / 



and 



Alu^j 2 -o; 2 )(4a; 2 -7 2 ) 
-2w 3 + (1 + 7 2 + 4(3 2 sin 2 §)w]£ 



(18) 



(19) 



where the function Q = Q(io, k) is given by: 

G = 16u; 4 -4(l+7 2 )w 2 +7 2 -4 / 3 2 (4w 2 -^ 2 )sin 2 fc. (20) 
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FIG. 4: (Color online) The dispersion coefficient P (top panel) and 
the nonlinearity coefficient Q (bottom panel) as functions of the fre- 
quency / (in GHz) in the LH band. Solid (blue) and dashed (red) 
lines correspond to the quasi-discrete and continuum approxima- 
tions, respectively. Note that there exist three characteristic frequen- 
cies, namely a frequency where the dispersion coefficient P van- 
ishes, and two frequencies corresponding to poles of the nonlinearity 
coefficient Q; on the contrary, in the continuum approximation, the 
coefficient Q possesses a sole pole. 



C. The quasi-discrete soliton solution 

As is well known (see, e.g., Ref. [ 46 ]), the NLS equation 
possesses soliton solutions, the type of which is governed by 
the signs of the dispersion and nonlinearity coefficients. In 
particular, if PQ > 0, the NLS model supports bright soli- 
ton solutions, while for PQ < it supports dark soliton 
solutions. In the case under consideration, the signs of the 
coefficients P and Q for Eq. 



17 1 depend on the frequency. 



In Fig. |4] we show this frequency dependence of P and Q 
in the LH frequency band (the figure shows P and Q both 
in the quasi-discrete and the continuum approximation, cor- 
responding to solid and dashed lines, respectively). First, it 
can readily be observed that, in the quasi-discrete approxima- 
tion, the dispersion coefficient (see top panel of Fig. |4]i satis- 
fies P < (P > 0) for / < 2.217 GHz (/ > 2.217 GHz) 
inside the LH frequency band. On the other hand, the nonlin- 
earity coefficient Q (see bottom panel of Fig. |4}, which has 
two poles (for / = 2.116 GHz and / = 2.199 GHz) inside 
the LH frequency band is Q > for / € [2.102,2.116) U 
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FIG. 5: (Color online) The signs of the product of the dispersion and 
nonlinearity coefficients, P and Q, for various subintervals of the LH 
frequency band within the framework of the quasi-discrete approxi- 
mation. The solid (blue) curve corresponds to the dispersion relation 
(see inset in Fig. [3}; the dashed (black) lines at / = 2.116 GHz and 
/ = 2.199 GHz indicate the location of the poles of the coefficient Q 
(see bottom panel of Fig. |4j, while the one at / = 2.217 GHz shows 
where the coefficient P vanishes (see top panel of Fig. |4). Notice 
that for PQ > (PQ < 0) the NLS model supports bright (dark) 
soliton solutions. 



(2.217,2.272] and Q < for / G (2.116,2.199) (in GHz). 
Thus, the product PQ takes the following signs: PQ > 
for / G (2.116,2.199) U (2.217,2.272], while PQ < for 
/ G [2.102,2.116) U (2.199,2.217) (in GHz); the above re- 
sults are summarized and demonstrated in Fig. [5] 

Here, it should also be noticed that in the continuum ap- 
proximation the nonlinearity coefficient Q has only one pole 
(at / = 2.21 GHz); as a result, this approximation estimates 
a broader interval wherein bright solitons can be formed: as 
may be seen in Fig. [4] (see also Appendix C), bright soli- 
tons can be formed in the intervals / G [1.976,2.203) U 
(2.21,2.272] (in GHz). While this pole arises close to the 
point where a pole emerges in the quasi-discrete approxima- 
tion, the other pole of the latter seems to be missed by the pre- 
sumably less accurate genuinely continuum approximation. 
In any case, the above investigation shows that bright soli- 
tons (corresponding to PQ > 0) are easier to be observed 
in the SRR-CPW structure: first, unlike dark solitons (cor- 
responding to PQ < 0), bright solitons are supported in a 
wide subinterval of frequencies inside the LH frequency band; 
on the other hand, they can be formed for carrier frequen- 
cies sufficiently far away from frequencies where the disper- 
sion coefficient P vanishes or the nonlinearity coefficient Q 
has resonances (see Fig. [4]). This is where we expect this 
type of approximations to be most relevant in describing the 
fully discrete dynamical system. Thus, below, we will con- 
fine our considerations to the case of bright soliton solutions 
of Eq. ( 17 1; a single soliton solution can be expressed in the 



following form (see, e.g., Ref. [46]): 

Vi(X,T) = r)sech[r](X - X (t))]exp[i(KX - CIT)}, (21) 

where -q is the soliton's amplitude (and inverse width), 
X (t) = X o (0) + VT is the soliton center, X o (0) is the ini- 



tial soliton location, V is the soliton velocity, K = V denotes 
the soliton's wavenumber, while 51 = (1/2)(J\~ 2 — rf) is the 
soliton's frequency. The above expression can be used to ap- 
proximate the unknown voltage V n (t) in Eq. ( 1 1 1, in terms of 
the original coordinates, as follows: 

V n {t) fa Vo sech[eri(n — cot — no)] cos(K n — Slot), (22) 

where no is the initial soliton position, while the soliton's am- 
plitude Vq, velocity cq, wavenumber Kq, and frequency SIq, 
are respectively given by: 



Vo 

CO 

Ko 



2er}y/2\PQ- 
v g + 2eKP, 
k + eK, 
uj + e(Kv g - 



2efiP). 



(23) 
(24) 
(25) 
(26) 



We conclude this Section by mentioning the following. The 
presented approximate soliton solutions in the LH frequency 
band have a unique feature as compared to ones that can be 
formed in the RH frequency band (i.e., for 3.395 GHz < 
/ < 4.734 GHz - see Fig. [3j. This can be understood by 
the fact that while in the RH regime the frequency increases 
with the wavenumber, in the LH regime the transmission line 
exhibits "anomalous" dispersion, and the frequency decreases 
with the wave number, as seen in Fig. [3] Thus, the solitons 
in the LH frequency band are actually backward waves, with 
group and phase velocities that are antiparallel to each other 
(see also discussion and relevant experimental observations in 
Refs. 0111321). 



m. NUMERICAL RESULTS 

A. Analysis of the "regular" SRR-CPW structure 

Let us now proceed to study numerically the evolution of 
the quasi-discrete solitons presented in the previous section 
in the framework of the fully discrete model of Eq. ( flT) . 
First, we will study the latter, taking parameter values as 
per Eq. ( 13 1 (recall that these were borrowed from Ref. 071 , 
which reported realization of this - as characterized here - 
"regular" SRR-CPW structure) and in the next Section, we 
will study experimentally relevant modifications in the SRR- 
CPW structure, so as to investigate how pertinent parameter 
changes affect the domains of existence of quasi-discrete soli- 
tons. 

In the results below, we have fixed the parameters related 
to the soliton's amplitude as r/ = 1 and e = 0.02, which cor- 
respond to an initial value of the voltage equal to Vo = 0.5 V 
(similar, and even smaller, values have also been used in ex- 
periments [29 -34 1). Furthermore, we have fixed the initial 
soliton position to X o (0) — 1/2, and we have varied the fre- 
quency / and the soliton wavenumber K (recall that the latter 
sets the initial soliton momentum). The chosen intervals of 
variation have been / G [2.11, 2.18] U [2.21, 2.25] (in GHz) 
and K G [0, tt]. Those values of / lie in the LH regime with 
both P < and Q < 0. A ring of N = 1001 cells (and 
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periodic boundary conditions) has been chosen for the simu- 
lations. 

In order to characterize the outcome of the simulations, and 
compare analytical and numerical results, we have defined 
two diagnostics. The first one is the evolution of the center 
(alias pseudo-center-of-mass), 



X{t) 



En nV n 

y v 2 



(27) 



and the second one is a measure of the width (alias pseudo- 
width), defined as 



W(t) 



(28) 



According to the results of the previous Section, analytical 
forms for these quantities can readily be found in the continu- 
ous setting [i.e., approximating the soliton as per Eq. (f2T|)]: 



X(t) 
W(t) 



X o (0) + c Q t 



2\/3e?7 ' 



2 + c at, 



(29) 
(30) 



Below we provide the outcome of some typical simulations 
(see Figs.|6]|9]l through density plots of V n , the spatial profile 
of V n at t = 2000, as well as the time evolution of the cen- 
ter of mass X(t) and the width parameter W(t). Generally, 
as we show below in more detail, the direct numerical inte- 
gration of Eq. ( [TT| ), with initial conditions borrowed from the 
analytical expression of Eq. (22 1 (for t = 0), have revealed 



the following: the quasi-discrete bright solitons exist, indeed, 
in the predicted frequency regions inside the LH frequency 
band; furthermore, their form, as well as the evolution of their 
center and width, can be well approximated by pertinent ana- 
lytical expressions provided above, especially in cases where 
the carrier frequency is chosen sufficiently different from cer- 
tain characteristic frequencies (i.e., where the dispersion co- 
efficient P vanishes or the nonlinearity coefficient Q has res- 
onances). 

Before proceeding with the description of our results, we 
should also note the following. Although most of our simula- 
tions were performed for relatively large normalized times - 
typically of order of t ~ 10 3 - given our time normalization, 
the physical unit time (set by the frequency f s h = 3.39 GHz) 
is very small, namely to = (27r/ s ? l )" 1 s» 50 picoseconds (see 
Sec. II.A). In fact, since all characteristic frequencies of the 
system (see Eq. ( 10 1) are in the microwave regime, all charac- 



teristic times are less than a nanosecond, rendering long sim- 
ulations extremely time-consuming. Nevertheless, in a par- 
ticular case where the condition above is fulfilled, we have 
performed a few extremely long simulations (with normal- 
ized time horizons of t = 10 7 , corresponding to a physical 
time of the order of a millisecond), finding that the agreement 
reported below is still upheld in these runs. This indicates that 
our predictions concerning soliton formation and robustness 
may be valid for experimentally relevant times. 

Let us expose our results starting with Fig. [6] which shows 
the case of a quasi-discrete soliton with carrier frequency 
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FIG. 6: (Color online) Top left: density plot of the time evolution of 
V n obtained numerically. The top right panel compares the analytical 
and numerical profiles of V n at t = 2000. The bottom panels show 
the time evolution of the center of mass (left) and the width diagnos- 
tic (right); in both cases, the solid line corresponds to the numerics 
and the dashed line to the analytical prediction. The parameters used 
are / = 2.15 GHz and K = 0. 



/ = 2.18 GHz and zero initial momentum, K = 0, which 
evolves as a stable object over long times. In this case, the 
agreement between analytical and numerical results pertain- 
ing to the soliton profile, but also to the evolution of the 
center of mass and width diagnostics, is very good. On the 
other hand, Fig. [7] shows the evolution of a soliton with also 
/ = 2.18 GHz but with nonzero soliton momentum, K = it 
(with similar conclusions), while the soliton of Fig. [8] corre- 
sponds to / = 2.11 GHz (for K — 0), which is close to 
the resonance of the nonlinearity coefficient Q. In this last 
case, it is clear that although quasi-discrete bright solitons ex- 
ist, the agreement between analytics and numerics becomes 
worse (especially as concerns the estimation of the soliton 
width parameter shown in the bottom right panel of Fig. IS). 
This can be attributed to the proximity to the resonance where 
we expect the conditions for the quasi-discrete approximation 
to be violated. Note that similar results have also been ob- 
tained for the upper (second) allowable regime inside the LH 
band [/ G (2.217, 2.272] (in GHZ)]. 

Finally, using parameter values corresponding to the case 
shown in Fig. |6](i.e., for / = 2.15 GHz and K = 0), we have 
performed a very long simulation, up to normalized times 
t = 10 7 (physical time ~ 1 ms) in order to check numerically 
the soliton robustness; details are shown in Fig. [9] As far as 
the width parameter W is concerned (bottom right panel of the 
figure), it is clear that - despite its fluctuations - it possesses 
an almost constant mean value which (although not captured 
precisely by the analytical approximation) indicates that the 
initial pulse does not spread out. Furthermore, the top pan- 
els of the figure - and particularly the snapshots of the pulse 
profile at t = 10 7 - clearly indicate that the soliton persists 
as a stable object up to this long simulation time (where the 
numerical integration was stopped). These results, as well as 
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FIG. 7: (Color online) Same as Fig. [6| but for / = 2.15 GHz and 

K = 7T. 
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FIG. 8: (Color online) Same as Fig. [6] but for / = 2.11 GHz and 

K = 0. 



the ones presented above, indicate that quasi-discrete solitons 
can be formed in SRR-based CPW nonlinear structures, and 
may propagate over experimentally relevant times. 



B. Modifications of the "regular" SRR-CPW structure 



In the previous section, we have studied in detail quasi- 



discrete solitons of Eq. (Hi that can be formed in the SRR- 



CPW structure for parameter values given in Eq. ( 13 I. Here, 
we will study different scenarios arising from modifications 
of either the CPW structure or the SRR geometry. As it is 
clear from Fig. [T] (and also the discussion in Sec. II.A and 
Appendix A), there is a considerable degree of flexibility as 
concerns the choice of the parameter values that are involved 
in the determination of the parameters 7, /1 and j3: indeed, one 
may consider different characteristic widths hi, h 2 of the di- 
electrics (and values of the relative dielectric constants e r i, e r2 



FIG. 9: (Color online) The soliton of Fig. [6j is evolved until t = 
2 x 10 5 . All the panels are similar to that of Fig. [6] except for the 
top right. In this panel, snapshots of the soliton at t = 2 x 10 6 and 
t — 10 7 are compared to the initial condition of the simulation in 
order to examine its robustness under a very lengthy time evolution. 



thereof), or different geometrical characteristics of the CPW 
structure or of the SRRs. 

Nevertheless, since the considered SRR-CPW structure has 
already been optimally realized in practice 11351 [37], one 
should consider changes that keep the basic characteristics of 
the configuration as close to its experimental realization as 
possible. We thus choose to keep the characteristics of the 
"main" dielectric substrate (characterized by the parameters 
hi and e r {), as well as the transverse width a of the CPW 
structure fixed. Furthermore, we also keep fixed the charac- 
teristics of the nonlinear dielectric (with parameters h 2 and 
e r 2), since - according to our considerations - they do not 
significantly affect the linear response of the system. On the 
other hand, we consider certain changes that would arise from 
a slightly different realization of the considered setup: in par- 
ticular, we will study the changes in the the width G of the 
slots for the CPW structure, and changes of the geometric 
characteristics of the SRRs, namely their radius r (for fixed 
width w and spacing d between the SRRs). As we will show 
below, these changes may improve the nonlinear SRR-CPW 
configuration, in the sense that they lead to an increase of the 
width of the LH frequency band, and also increase the "cen- 
tral" frequency band (i.e., in between the resonances of co- 
efficient Q) where bright solitons can be formed (PQ > 0). 
Other changes in the parameter values have also been stud- 
ied, but the results will not be exposed here, as they lead to 
results qualitatively similar to the ones presented above, or 
even worse (i.e., they lead to decrease of the widths of the 
LH frequency band and/or the domains of existence of bright 
solitons). 

A relevant study along the above lines has led to the fol- 
lowing results. First, for a fixed value of the distance between 
the slots, D = 7 mm, we have changed the value of the slot 
width G in the interval 0.7 mm — 1.5 mm (recall that originally 
G = 1.35 mm). Such changes result in different values of Kj 
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FIG. 10: (Color online) Top panel: Same as Fig. [5] but for a slot 
width G = 1 mm. Middle and bottom panels: Same as Fig. [4] (for 
G = 1 mm). In this case, the values of parameters 7, fj, and j3 are 
given in Eq. {3T). 



(j = 0, 1, 2) in Eq. ( A4i and, accordingly, to different values 
of the effective circuit model parameters Cr, Cl, Ll, C s and 
L s , which, in turn, provide different values for the parameters 
7, /i and j3 [cf. Eq. (10l]. As G is decreased (increased) in 
the aforementioned interval, the parameters 7 and /i are also 
decreased (increased), while the parameter f3 is increased (de- 
creased). These changes result in a decrease (increase) of the 
LH frequency, but also lead to an increase (decrease) of the 
"central" frequency band where bright solitons can be formed. 
As an example, in the top panel of Fig. 



10 we show the LH 



frequency band, as well as the regimes for soliton formation, 
in a case corresponding to G = 1 mm, for which the charac- 
teristic parameters of the configuration take the values: 




FIG. 11: (Color online) Same as Fig. [9] but for parameter values 
given by Eq. \3l\ . Notice that the soliton is evolved until t = 2 x 10 4 
(top left panel), but it persists for much longer times, namely up to 
t = 10 7 (see snapshots in the top right panels). Once again, the 
width parameter W (bottom right panel) possesses an almost con- 
stant value, which indicates the soliton-like nature of the initial pulse. 



As seen in this figure, the LH band extends from 2.095 GHz to 
2.245 GHz (i.e., it is decreased by ps 11%), while the central 
frequency region, where bright solitons exist, extends from 
2.103 GHz to 2.201 GHz (i.e., it is increased by « 15%). No- 
tice that the functional form of the parameters P and Q shown 
in the middle and bottom panels of Fig. [10] respectively, are 
similar to the ones in Fig. [4] One should notice, however, that 
the parameter Q exhibits an almost flat profile in the central 
frequency band with PQ > 0. 



Numerical simulations shown in Fig. 1 1 for bright soliton 
propagation at the frequency / = 2.15 GHz lead to results 
qualitatively similar to the ones presented above (cf. Figs. [6] 
and [9]). In particular, once again, the soliton is quite robust 
up to long times (see the snapshots corresponding to t = 5 x 
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while the 
has an 



7 = 0.61, n = 0.54, 



P = 0.48. 



(31) 



10 6 and t = 10 7 in the top right panel of Fig 
width parameter W (bottom right panel of Fig 
almost constant mean value (close to the value corresponding 
to the analytical estimate) indicates that the initial pulse does 
not spread out, thus featuring genuine soliton characteristics. 

Next, keeping fixed the parameter values of the CPW struc- 
ture as in Sec. II. A (in this case G = 1.35 mm), we have stud- 
ied changes in the SRR geometry, namely the radius r and 
width w of the SRRs (for fixed spacing between the SRRs, 
d = 0.2 mm). First, in the case with fixed w = 0.6 mm, 
the increase (decrease) of the SRR radius r - in the interval 
1.8 mm to 1.6 mm - results in a decrease (increase) of the 
values of parameters 7 and \i, while the parameter (3 takes the 
approximately constant value of j3 w 0.48. Accordingly, it is 
found that the width of the LH frequency band is increased 
(decreased), while the central frequency regime where bright 
solitons can be formed is increased (decreased) as well. As 
an example, for r = 2.6 mm (recall that r = 2.4 mm for the 
"regular" SRR-CPW structure), the LH frequency band ex- 
tends from 2.106 GHz up to 2.343 GHz, i.e., it has increased 
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~ 39%. Furthermore, in the same case, the central frequency 
region where bright solitons can be formed is also increased 
by s» 29%, with the frequency dependence of coefficient Q 
being similar to the one shown in the bottom panel of Fig. 10 
In this case, the parameters 7, /i and j3 take the following val 
ues: 



7 = 0.69, n = 0.57, (3 = 0.48. 



(32) 



Numerical simulations for bright soliton propagation in this 
regime (e.g., we have performed long evolution runs for the 
frequency / = 2.181 GHz) lead to results qualitatively similar 



to the ones presented above (cf. Figs.kM^and 1 1 



IV. DISCUSSION AND CONCLUSIONS 

We have studied analytically and numerically quasi- 
discrete microwave bright solitons that can be formed in the 
left-handed frequency band of a split ring resonator (SRR) 
based coplanar waveguide (CPW) structure. We have used the 
nonlinear transmission line analogue of this structure to de- 
rive a nonlinear lattice equation governing the voltage across 
the fundamental (unit cell) element of the system. This lat- 
tice equation was then treated analytically, by means of the 
quasi-discrete approximation. The latter is a variant of the 
multi-scale perturbation method, which takes into regard the 
discreteness of the system by considering the carrier (enve- 
lope) of the wave as a discrete (continuum) object. This ap- 
proach allowed us to derive, in the small-amplitude approx- 
imation and for certain space- and time-scales, a nonlinear 
Schrodinger (NLS) model for the unknown voltage envelope 
function. The NLS model was then used to predict formation 
of bright solitons in certain sub-regions of the left-handed fre- 
quency band. Importantly, the conditions for soliton forma- 
tion were found to depend on the discreteness of the system: 
in fact, if the continuum - instead of the quasi-discrete - ap- 
proximation was used, then the allowable (for soliton forma- 
tion) frequency bands would be significantly modified. Fur- 
thermore, the quasi-discrete approximation predicts effects, 
such as the appearance of resonance frequencies in the coeffi- 
cients of the effective NLS model, which suggest optimum op- 
erating frequency bands for the observation of quasi-discrete 
solitons. Generalizing these observations, one should ex- 
pect that the adopted analytical approach, based on the quasi- 
discrete approximation, should provide more reliable results 
concerning conditions for soliton formation in left-handed 
structures. 

Numerical simulations performed in the framework of the 
nonlinear lattice equation, with initial conditions borrowed 
from the effective NLS equation, revealed that the (bright) 
solitons may indeed be formed in the frequency sub-intervals 
predicted in the analytics. Furthermore, it was shown that if 
the carrier frequency is chosen to be sufficiently far from char- 
acteristic frequencies (where the dispersion and nonlinearity 
coefficients of the effective NLS model vanish or have reso- 
nance poles), the numerically found soliton profile and char- 
acteristics (center of mass and width) are in very good agree- 
ment with the analytical predictions. Although the numeri- 



cal simulations were performed for relatively large computa- 
tion times, the corresponding physical times were small due 
to the fact that all the characteristic frequencies of the system 
were in the microwave regime: under any choice, the physi- 
cal time unit is less than a nanosecond and, thus, simulations 
corresponding to a few seconds would correspond to compu- 
tation times of order of 10 10 . This difficulty, particular to the 
microwave structure under consideration, did not allow us to 
be definitively conclusive as concerns the robustness of the 
solitons for realistic experimental times. Nevertheless, very 
long simulations (corresponding to computing times ~ 10 7 or 
~ 1 millisecond in physical units) have shown that the pre- 
dicted quasi-discrete microwave solitons are quite robust and 
do not appear to be modified during their propagation. This, in 
turn, allows us to conjecture that these nonlinear wave struc- 
tures have a good chance to be observed experimentally in the 
near future. 

We have also presented a study showing how certain phys- 
ical parameters of the configuration may affect the results. In 
particular, we have focused on two cases, namely the effects 
of decrease of the width of the slots in the CPW structure, 
and of the increase of the SRR radius. These choices led to 
modified values of the parameters of the model (/?, 7, and /i), 
such that the width of the left-handed frequency band or/and 
the domain of existence of quasi-discrete bright solitons are 
increased. This way, we have also proposed certain experi- 
mentally relevant changes in the considered configuration, to 
facilitate observation of solitons in future experiments. 

There are many interesting directions for future studies. In 
that regard, first, we should mention that in the analysis of the 
considered SRR-CPW structure we have actually excluded the 
SRR from nonlinearity. It would be interesting to see what 
happens if SRR is also subject to nonlinear modulation: this 
is certainly a very challenging direction, in terms of the de- 
velopment of both analytical and numerical techniques, due 
to the form of the pertinent models, which are coupled non- 
linear lattice equations [see, e.g., Eqs. ^-{jj}]. Such studies 
would lead to potentially interesting and relevant results con- 
cerning nonlinear wave propagation in such settings, as well 
as the design of nonlinear left-handed transmission lines (and 
related structures). 

On the other hand, we note that following the analytical 
procedure adopted in this work, it would be possible to ana- 
lyze soliton formation in relevant left-handed structures. Fur- 
thermore, employing this analytical approach, one could - in 
principle - derive self-consistently [at 0(e 4 ) in the perturba- 
tion scheme] a generalized NLS model, incorporating higher- 
order dispersive and nonlinear effects (see, e.g., Ref. [24]). 
Such a higher-order NLS model could also predict soliton 
propagation but, in this case, the soliton characteristics would 
be modified by the presence of the higher-order effects. Thus, 
an interesting challenge would be the derivation of such a 
higher-order NLS model, the study of its soliton solutions and 
a comparison of such findings to direct simulations. 

Finally, while the theory and simulations presented in this 
work assume ideal components, for which excellent perfor- 
mance can be observed, the reality to be met in possible ex- 
periments might deteriorate the performance by unavoidable 
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dispersion in component parameters. The effect of disorder is 
expected to be quite significant in nonlinear settings: in fact, 
its paramount importance has been demonstrated even in the 
simpler case of linear resonant systems (see, e.g., Ref. fl7l ) 
and, thus, the presence of disorder may drastically affect the 
results in nonlinear structures. In that regard, a study of how 
much the considered system (or other relevant ones) is sensi- 
tive to disorder is certainly a relevant and important direction 
for future studies. 
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In the above expressions, F « 0.54 is the fractional area of 
the slots occupied by the rings ll37ll . while Lsrr an d Csrr 
denote, respectively, the inductance and capacitance of the 
SRRs, and are given by: 
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where /io is the permeability of vacuum, r, d and w respec- 
tively denote the radius of the internal ring, the distance be- 
tween the internal and external rings, and the width of the 
rings (see Fig.[TJ), while the arguments of the elliptic integrals 
K and K' are now given by: 



Appendix A: The physical parameters of the SRR-CPW 
structure 

The elements Cr, Lr and Ll associated with the CPW 
structure are connected with the physical parameters of the 
system by means of the following equations (see, e.g., 
Ref. HI): 
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In the above expressions, eo is the dielectric constant of vac- 
uum, K and K' are the complete elliptic integral of the first 
kind and its complementary function, respectively 11441 . the 
arguments of these functions are given by: 
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and, finally, L p and Lcpw = 4e (307r) 2 K' '(kq) / 'K(ko) de- 
note, respectively, the inductance of the shunt strip and the ef- 
fective inductance of the CPW structure (note that the former 
takes the approximate value L p ss 0.36 nH [37 1). As far as 
the values of the SRR parameters, L 8 and C s , are concerned, 
they are given by the following expressions |[37l : 
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- i 7 ' 2 L^p W /2LpLsRR ' 

1 + F 2 L CPW /2L p Lsnn\ 



l + L 



CPW 



(A5) 



(A6) 



sinh 



k-3 



d + 2w' 



sinh 



Tr(d+2w) 

IF, 



(A9) 



Appendix B: The perturbation scheme 



As mentioned in Sec. II. B, the quasi-continuum approxi- 
mation is a variant of the method of multiple scales B31 . We 
thus introduce, at first, a set of new independent temporal vari- 
ables, t n — e n t (n = 0, 1, 2, • • • ), and acordingly expand the 
derivative operator d t as d t — dt + edt 1 + ■ ■ ■ ■ Next, we seek 
solutions of Eq. (JTTJ in the form: 



V n = eu ln (t n )e i0n + e 2 U2n(t n )e rM " +■■■ + 



211 „ 



(Bl) 



where 8 n = ujta — kn. Then, we substitute Eq. ( |B1| > into 
Eq. ( [TT] i and employ a continuum approximation for the un- 
known envelope functions u n , namely u n u(x), where 
x = na and a being the lattice spacing (the latter parame- 
ter does not appear in the results below, as one may readily 
rescale x as x/a). Furthermore, similarly to the introduction 
of the temporal variables, we introduce the set of the spatial 
variables x n = e n x (and, thus, d x = d Xo + ed Xl +•••)■ To 
this end, equating coefficients of like powers of e, we obtain 
the following (first three) perturbation equations: 



0(e) 
0(e 2 ) 
0(e 3 ) 



L ux = 0, 
L u 2 + L\u\ 
Liu 2 + L 2 u 2 



N ul, 
Nohnuz] 



(B2) 
(B3) 
(B4) 
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where the operators are given by 



Ln = 



(f_ 

04 



l +7 2 + 4/Vsin 2 - 



4/3 V sin 2 -+ 7 2 , 



(B5) 



U = 



+ 2(1 + 7^+ 4/3 2 /i 2 sin^ - 



if 



2i/3 2 sinfc 



dtldxi ' 



G 



k 

2 ) dtodh 

(B6) 

/>■ \ d 2 



^fe +2 ( 1 + 72+4/3Vsin2 2/'^ 



iV 



4i/3 sinfc 

d 4 



2i/3 2 ^ 2 sin k 



d 2 



dxi ' 



(B7) 



(B8) 



It is clear that the first-order pertubation solution to Eq. ( B3 i 
reads: 

ui = Vi (xx, X2, ■ ■ ■ , ti , *2, • • • ) exp(i#) + c.c. , (B9) 
where V\ is an unknown complex function, 8 — ujto — kxo, 



while u> and k satisfy the dispersion relation of Eq. ( 15 1. Next 



substituting Eq. ( B9 1 into Eq. ( B4 1, we obtain the following 



results: first, the non-secularity condition: 



dVi 



f3 2 smk(u 2 -^ 2 ) 



2u 3 - (1 + 7 2 +4/3 2 sin 2 §)w 



dVi 
dx\ 



0, (BIO) 



suggests that V\ = V\(X, X2, • • • , £2, • • • ), where X — X\ — 
Vgti and v g is the group velocity of Eq. ( [T6| , now consistently 
determined; second, we obtain a uniformly valid solution for 
the second-order perturbation equation, in the form: 



"2 



4w 2 ( 7 2 -4w 2 ) 



exp(i29) + c.c. 



(Bll) 



where Q{ui,k) is given by Eq. (20i and the dependence of u-i 



on higher-order scales has been omitted. To this end, substi- 



tuting Eqs. (B9 i-(Bl 1 1 into Eq. (|B4[), and using the variables 



X = x\ — Vgti = e(n — Vgt) and T = <2 = e t, we derive 



from the non-secularity condition at 0(e 3 ) the NLS Eq. ( 17 1. 



Appendix C: The continuum approximation 

In the continuum limit (for k —> 0), the dispersion relation 
[cf. Eq. ( [7~5] >] is reduced to the form, 

- (1 + 7 2 + f3 2 k 2 )oj 2 + f3 2 n 2 k 2 + 7 2 = 0. (CI) 

The group velocity, v g = dui /dk, is now given by 

_ /? 2 fc(^-M 2 ) rr9 , 
9_ 2^-(l+7 2 + /3 2 fc 2 V { ' 

Finally, the expressions for the dispersion and nonlinearity co- 
efficients P and Q in the continuum approximation read: 



P 



and 



(1 + 7 2 + fi 2 k 2 - 6lu 2 )v 2 
-4w 3 + 2(1 + 7 2 + P 2 k 2 )uj 

4p 2 ukv g + P 2 {lu 2 - n 2 ) 
-4w 3 + 2(1 + 7 2 + ^ 2 A; 2 )a;' 



4w 4 ( 7 2 -uj 2 )(4uj 2 - 7 2 ) 



(C3) 



(C4) 



^ [-2uj 3 + (1 +7 2 +f3 2 k 2 )w]g' 
where the function Q is now given by: 

g = 16w 4 - 4(1 + 7 2 )w 2 + 7 2 - 4^ 2 (4cj 2 - fi 2 )k 2 . (C5) 
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